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2 Abstract 

i— i Current Higgs data at the Large Hadron Collider is compatible with a SM signal at the 2a level, 

but the central value of the signal strength in the diphoton channel is enhanced with respect to 
^s! the SM expectation. If the enhancement resides in the diphoton partial decay width, the data 

could be accommodated in the Minimally Supersymmetric Standard Model (MSSM) with highly 
mixed light staus. We revisit the issue of vacuum instability induced by large mixing in the stau 
sector, including effects of a radiatively-corrected tau Yukawa coupling. Further, we emphasize 
the importance of taking into account the tan/3 dependence in the stability bound. While the 
metastability of the Universe constrains the possible enhancement in the Higgs to diphoton decay 
width in the light stau scenario, an increase of the order of 50% can be achieved in the region of 
large tan/3. Larger enhancements may be obtained, but would require values of tan/3 associated 
with non-perturbative values of the tau Yukawa coupling at scales below the GUT scale, thereby 
implying the presence of new physics beyond the MSSM. 
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I. INTRODUCTION 



The discovery of a new boson at the CERN Large Hadron Collider (LHC) is an extraor- 
dinary achievement for high energy physics [TJ [2]. Preliminary studies indicate that the 
properties of the new boson are generally consistent with that of a Standard Model (SM) 
Higgs boson [3J. It is of course crucial to eventually measure the spin, CP, and electroweak 
quantum numbers of this new particle. In current data the significance of the discovery 
is mainly driven by the observation of an excess of events in the diphoton and four-lepton 
channels, that are compatible at the 2 a level with a SM-like Higgs boson with a mass at 
around 125 GeV. 

Additionally, there are indications that the new particle decays into WW with SM-like 
rates. Searches are also being performed in the bb and rr channels. The decay rates in these 
channels also show consistency with SM expectations at the 1 a level |4j,|5j. Although the 
diphoton signal is compatible with a SM Higgs signal at the 2 a level, the central value of 
the signal strength is observed to be approximately 2 and 1.5 times the SM prediction by 
the ATLAS and CMS collaborations, respectively [6]. Taken at face value, the data seems 
to suggest that the enhancement in the diphoton rate arises from an enhanced partial decay 
width of the Higgs to diphotons, while the production cross-sections seem to indicate that 
the couplings of the new resonance to vector bosons and top quarks are similar to those 
predicted for the SM Higgs. At present, the measurements are statistically limited, and 
more data is necessary to reach conclusive results. 

The deviation in the diphoton rate from the SM prediction, if confirmed, would be a 
clear indication of physics beyond the SM. In the MSSM, the best motivated candidates 
that can significantly modify the Higgs to diphoton decay width are sleptons, in particular 
staus [3 E]. A significant enhancement requires that one of the staus is light, of the order 
of 100 GeV, and further that the mixing between the left- and right-handed staus is large. 

More generally, an increase in the Higgs to diphoton partial width without impacting the 
production rate requires new color neutral, electrically charged particles with masses of the 
order of 100 GeV and significant couplings to the Higgs boson [7]-[TT]. Such new particles 
inevitably modify the Higgs potential through quantum corrections. If the new charged 
particles are fermions, they may drive the Higgs quartic coupling negative via renormal- 
ization group evolution and thus destabilize the Higgs vacuum at scales of the order of a 
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few TeV [10]. For new charged scalars, an enhanced diphoton width requires a large cubic 
coupling between the Higgs and a pair of the new scalars, which in turn could induce a 
new charge-breaking vacuum [5J. It is therefore important to study the extent to which 
the diphoton partial width could be enhanced via light charged particles without inducing 
vacuum instability at energies well below the GUT scale. 

In the MSSM, constraints from charge-breaking minima induced by large stau mixing 
were studied in Ref. [12], and later refined in Ref. [13]. More recently, Refs. [HI [15] consid- 
ered vacuum stability issues in the heavily mixed, light stau scenario with enhanced Higgs 
to diphoton partial width based on the results in Ref. [13J. In this work we revisit the 
vacuum stability issue and improve upon the work in Ref. [13]. In particular, we emphasize 
the importance of retaining the dependence on the ratio of the Higgs vacuum expectation 
values, tan/3, for fixed values of the stau mixing mass parameter, /xtan/3, when deriving 
vacuum stability bounds. We also include important tan /3-enhanced effects from radiatively 
corrected tau Yukawa couplings. Our study shows that a diphoton partial width enhance- 
ment of 0(50%) above the SM expectation can be compatible with the metastability of 
the electroweak-breaking vacuum for sufficiently large tan /3 in the MSSM. An enhancement 
beyond (9(50%) requires such large values of tan/3 that the r Yukawa coupling becomes 
non-perturbative below the GUT scale, which in turn would imply new physics beyond the 
MSSM at scales below the corresponding Landau pole. 

This work is organized as follows. In Sec. |TT] we analyze the dependence of the vacuum 
stability conditions on the parameters of the model, such as Yukawa couplings, A T , and tua, 
and compare with previous results. We further make connection with the enhancement of 



the Higgs to diphoton decay width with these same parameters. In Sec. II A we investigate 
the radiative corrections to the tau and bottom Yukawa couplings and their impact on the 



vacuum stability conditions. Sec. II B discusses the effects of A T and on the total width 
of the Higgs decay, which enters into the diphoton decay branching fraction. Numerical 



results, which follow from these considerations, are presented in Sec. II C We then discuss 



some of the possible constraints on the large tan/3 region in Sec. |III[ In Sec. [TV] we present 
our conclusions. 
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II. STAUS, HIGGS DIPHOTON WIDTH AND VACUUM STABILITY 



To study the metastability of the electroweak vacuum in the MSSM, in the presence of 
light staus with large mixing, it is instructive to first write down the scalar potential for 
the neutral component of the up-type Higgs, h u , the left-handed stau, fi, and the right- 
handed stau, tr. Neglecting the down-type Higgs is a very good approximation since for 
large values of tan /3 - as required to achieve large mixing in the stau sector- and sizable 
rriA, the vacuum expectation value (VEV) remains small. However, for completeness, 



when presenting most of our numerical results in Section II C we will use the full potential 
including both the up- and the down-type Higgs bosons. 

Following Ref. [32], and normalizing all fields as complex scalar fields, the scalar potential 
can be written as 

v = i/i^-^r i f*i 2 +|(if L r+i^r) 2 +|(if L | 2 -2i^i 2 -i/ iu | 2 ) 2 

o o 
i 2 it, |2 i 2 i~ |2 i 2 I "* 1 2 i 9l~^~92 r it |4 /i\ 

+m H Jh u \ +m L& \T L \ +m E3 \r R \ H — H \h u \ , (1) 

where /i is the Higgsino supersymmetric mass parameter and y T is the tau Yukawa coupling 
appearing in the MSSM superpotential. #2 and g\ are the gauge couplings for SU(2)l and 
U (l)y, respectively. In addition, m 2 H , m 2 L3 , and m E are the soft-breaking masses for the up- 
type Higgs, the left-handed third generation sleptons, and the right-handed third generation 
sleptons. The last term, proportional to 5h in Eq. ([I]), represents the leading contribution to 
the full one-loop effective potential, arising from the top and stop loops. This contribution 
depends on the average stop mass, m^, and on the stop mixing parameter, X t = A t — fi cot /3, 
and is of the order of unity for a 125 GeV Higgs boson [T3] , 
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where y t ~ \f2m t /v with v ~ 246 GeV, and m t is the weak scale running top-quark mass. 

The source of vacuum instability in Eq. ([T]) is clear: the term coupling the Higgs to the two 
staus, whose coefficient is proportional to fiy T , has a negative sign which tends to destabilize 
the vacuum for positive values of the fields. As first studied in Ref. [12], when this cubic 
coupling becomes too large, a charge breaking vacuum deeper than the electroweak breaking 
vacuum may exist. Moreover, after the Higgs acquires a VEV, h u (d) — > (tVd) + ^«(d)))/"\/2, 
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this cubic coupling, ny T /y/2, also contributes to the off-diagonal entry in the stau mass- 
squared matrix. Including tan /3-suppressed terms that were not included in the potential 
in Eq. Q, the stau mass matrix is given by 

M 2 = ( m 2 L3 +m 2 T + D L ^(A r cos/?-//sin/3)\ 
y 5§ (A T cos/3-/isin/3) m 2 E3 +ml + D R J 

where A T is the soft-breaking A-term for staus, and Dl r are the D-term contributions to 
the slepton masses. Therefore the coefficient that triggers the vacuum instability is also 
crucial in determining the mixing in the stau sector and hence the possible enhancement in 
the Higgs to diphoton width. 

A charge breaking minimum is not necessarily a problem if the ordinary electroweak 
breaking vacuum is metastable with a lifetime longer than the age of the Universe. The 
lifetime of a metastable vacuum is usually computed using semiclassical techniques. The 
probability of decaying into the true vacuum per unit spacetime volume is given by [16] 

£ = A e -*», (4) 

where A ~ (100 GeV) 4 is a dimensionful parameter, expected to be roughly the fourth 
power of the electroweak scale. The Se is the Euclidean action evaluated on the "bounce" 
solution that interpolates between the metastable vacuum and the other side of the barrier. 
The volume is given by V = R 4 , with R being the characteristic size of the bounce. Asking 
for the lifetime of the metastable vacuum to be longer than the present age of the Universe 
is equivalent to requiring T/V to be smaller than the fourth power of the Hubble constant, 
H ~ 1.5 x 10~ 42 GeV. This then implies that the vacuum is metastable if [T7J [18] 

Se > 400 . (5) 

In Ref. [T3] a numerical study based on Eq. (jlj in the configuration space of the three 
fields h u , f l and tr, found that the metastability condition is mainly sensitive to //tan/3, 
rriL 3 i an d TnE 3 , while the dependence on tan/3 for fixed /itan/3 and on 5h is small. The 
resulting vacuum metastability condition in Ref. [13] was summarized as follows 1 

\fi tan fi\ < 76.9y/m L3 m E3 + 38.7(m L3 +m E A- 1.04 x 10 4 GeV . (6) 



1 We recently learned that Eq. |6| is being revised [TH] 



The small dependence of the metastability condition on Sh, as claimed in Ref. [13], is under- 
standable since it only affects the Higgs quartic term. In particular, the bound on /ztan/3 
changes by only ~ 10% when Sh is varied from to 1. In our study, we will not investigate 
this dependence further and, unless otherwise stated, fix Sh at 0.9 since the mass of the 
Higgs has been measured to be approximately 125 GeV. 

The dependence on tan /3 can be understood starting from the tree-level relation for the 
tau lepton mass and the tau Yukawa coupling, 

y T ^V2 — 7 

v cos p 

When tan/3 > 1, sin/3 « 1 - l/(2tan 2 /3), which very quickly approaches unity, hence we 
can approximate sin/3 ~ 1. Then Eq. ([7]) implies that for tan/3 ^> 1, the dependence of y T 
on tan /3 is, to a very good approximation 

(8) 

v smp 100 

As mentioned earlier, the coefficient of the destabilizing cubic term in the scalar potential, 
huTL^R, is proportional to fiy T , and therefore to /itan/3. However, there is also a term 
which depends only on tan 2 /3 which arises from the stabilizing quartic term, |tl77j| 2 , whose 
coefficient is y 2 T . Therefore, for fixed /ztan/3, or equivalently fixed /j,y T , there is a residual 
dependence on tan /3 coming from this stabilizing term. 

Using the constraint in Eq. ^ from Ref. [13], that directly relates the vacuum metasta- 
bility condition to the stau mixing in the MSSM, the author of Ref. [15] derived a maximum 
allowed enhancement of the Higgs to diphoton width of about 25% of the SM value, for a 
light stau mass heavier than about 100 GeV. Our work, however, shows that the bound on 
/itan/3 is about 15% larger than the one displayed in Eq. ([6]). Moreover, we will show that 
the residual dependence on tan/3, for a fixed /itan/3, can have a significant impact on the 
vacuum metastability requirement. This is to be expected since if /itan/3 (or equivalently 
fiy T ) is held constant, then going to larger values of tan /3 increases the coefficient of the sta- 
bilizing quartic term, Ir^f^l 2 , and alleviates the vacuum stability constraint. This implies 
that larger values of tan /3 will lead to a further relaxation of the bound on //tan /3, allowing 
for larger enhancement of the Higgs to diphoton rate from light stau loops. As we will show 
below, within the MSSM the magnitude of the possible enhancement is thus constrained by 
how large the value of tan /3 (or equivalently of y T at low energies) can be, without implying 
that the tau Yukawa coupling develops a Landau Pole at energies below the GUT scale. 



We shall make comparisons with previous works, such as Ref. [13], using the improved 
tree-level approximation in Eq. (fiT). However, our final results are obtained using the full 
one-loop effective potential, including also terms involving either hd or suppressed by tan /3 
which were neglected in Eq. ([!]). Including such terms, in the large tan/3 limit, we obtain 
additional contributions to the scalar potential: 



AV ~ m\\h d \ 2 - ^n{hdh u + h.c.) + ^^\hu\ 2 + {y T A T h d f L f* R + h.c.) 

+ |2/t| 2 N 2 (|^l| 2 + \r R \ 2 ) + D-terms. (9) 

Comparing this with Eq. we see that and A T could also impact vacuum stability. 
First we note that for positive values of A T and //, in the charge breaking minima, the field hd 
tends to acquire values which are opposite in sign to the h u values. This means that positive, 
non-negligible values of A T contribute constructively to the destabilizing trilinear term in the 
scalar potential, thereby tightening the bound on // tan /3 from the metastability condition. 
However, the m 2 A terms still give a positive contribution to the scalar potential. Therefore, 
the A T effects are suppressed for large values of tua, which then tend to reduce the values 
of h d associated with the effective potential minima. Hence we see that additional charge 
breaking minima may be induced for small values of and large values of A T . This in 
turn implies that non-zero A T with \iA T > 0, can further suppress the possible enhancement 
of the diphoton width coming from light staus, depending on the value of tua- 

We are interested in regions of parameter space where A T is smaller than or of the order 
of 1 TeV and //tan/3 ~ 0(30) TeV. In this region, A T does not directly play an important 
role in the value of the stau mass, as can be seen from Eq. (J3|. However, for a given set 
of parameters, a positive value of \xA r lowers the minimum value of // tan /3 allowed by the 
metastability condition, suppressing the mixing effect in the stau sector. This then increases 
the lightest stau mass compatible with vacuum stability and reduces the enhancement in 
the diphoton partial decay width. On the other hand, positive values of \lA t reduce the 
total width of the lightest CP-even Higgs via mixing in the CP-even Higgs sector |7J [5] . It 
turns out that the reduction in the enhancement of the diphoton partial width is largely 
compensated by the decrease in the total decay width, leaving a diphoton branching fraction 
that is not severely affected by A T . The effect on the total width of the Higgs due to mixing 



effects in the Higgs sector will be discussed in detail in Sec. II B 
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A. Radiative Corrections to the Tau and Bottom Yukawas 



As stated in Ref . [13] , the reason the stabilizing effect from tan (3 is naturally suppressed 
is because the tree-level tau Yukawa coupling, y T , is proportional to tan/3 with a small pro- 
portionality constant ~ 1/100, Eq. rtsl) , and only very large values of tan/3 would overcome 
that suppression. The proportionality constant in Eq. ([7| is modified at one-loop level to 
be [20] 

- /o mT tan/3 

Vr ~ V t;cos/3(l + A T ) * 100(1 + A T ) ' UUj 

where A T arises dominantly from a stau-neutralino and a sneutrino-chargino loop, and may 
become of (9(0.1) for sufficiently large values of tan/3. An approximate expression for A r 



is given in the Appendix in Eq. (22). The full analytic expression for A r can be found in 
Ref. |20j and has been implemented in a new version of the code CPsuperH [22], which we 
use in our study. 



The bottom Yukawa coupling receives similar modifications, 

Vb 



V %cos/3(l + A 6 ) ' 1 J 

where A& is dominated by contributions from sbottom-gluino and stop-chargino loops. An 



approximate expression for A& is also give in the Appendix. From Eq. (23) we see that, 
for stop masses of the order of 1 TeV and sizable trilinear terms A t that are needed to 
accommodate a 125 GeV Higgs, the stop-chargino loop contribution becomes sizable. So, 
while the sbottom-gluino loop contribution could be suppressed by a large mass splitting 
between the sbottom and the gluino, the stop-chargino loop contribution in A& is always 
sizable in the scenario we consider 2 . 

The A T corrections are smaller in magnitude compared to the A& corrections because they 
are suppressed by electroweak gauge couplings. Additionally, A r tends to be dominated by 
loops which include electroweak gauginos and acquires a sign opposite to that of fiM 2 . On 
the other hand, for squark and gluino masses that are of the same order, the sbottom- 
gluino contribution to A& becomes the dominant one and A& acquires the same sign as 
yuM 3 . Furthermore, the stop-chargino contribution is proportional to fiA t and adds to the 

2 Alternatively, one could have stop and sbottom masses of the order of a few tens of TeV to obtain a 125 
GeV Higgs without a significant A t . In that case, both the gluino-sbottom as well as the stop-chargino 
loop would be suppressed and A& would be small. 
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effect of the sbottom-gluino loop if \xA t has the same sign as jjM^. In the following we 
shall consider values of \iMi > for all three gaugino masses, i = 1,2,3, and fiA t > 0, 
which then lead to positive values of Af, and negative values of A T . Note that this choice of 
signs improves the agreement between the theoretical prediction of the anomalous magnetic 
moment of the muon and its measured value [23] , [21] and also helps in weakening the bounds 
on the value of tan coming from the requirement of keeping the bottom Yukawa coupling 
perturbative until high energy scales. Moreover, as we shall discuss in more detail below, 
positive values of \xA t are helpful in avoiding the constraints coming from the B s — > fifi rare 
decay measurement 



From Eq. (10) it is clear that one can define an effective tan/3 in the tau sector 3 by 



n tan/3 

tan/3 reff = - - - , (12) 

which simplifies the relation between y T and m T . Looking at the stau mass matrix, given in 
Eq. (|3]), we see that if we further define an effective A T by 

A "« = TT^- (13) 

then the stau mass-squared matrix can be re-written as 

_ i m 2 L3 + m 2 T + D L m T (A TeS - /itan/3 reff ) ^ 
m T (A Tefi - /xtan/3 reff ) m| 3 + m 2 . + D R 

Since the relation between m T and y T , as well as the stau mass-squared matrix, retain their 



tree-level form when using the effective tan /3 defined in Eq. (12), we will find it convenient 
to express our results in terms of this effective quantity. 

In Fig. [l]we plot A T and y T as functions of tan/3 for ttil 3 — mE 3 = 250 GeV, fi = 520 
GeV, M 2 = 400 GeV and A T = 1 TeV. The shaded region in the figures corresponds to 
a stau mass that is below the LEP limit of about 90 GeV [26]. From the left panel we 
see that A T is negative and of the order of 10 - 15% for tan j3 = 40-60. The right panel 
of Fig. [T] shows the effect of A r in increasing the value of y T for any given value of tan/3. 
For example, for tan /3 = 50 the corresponding y T value, without the inclusion of the A r 
effect, can be read-off using the line labeled "Tree" to be ~ 0.51. Including A T increases 



3 Notice that one could define a different effective tan /? in the bottom sector. 
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FIG. 1: A r versus tan/3 (left panel) and y T versus tan/3 (right panel) for m£ 3 = rriE 3 = 250 GeV, 
/i = 520 GeV and A T = tua = 1 TeV. In the right panel, the red (solid) line uses the one-loop 
relation, including A T , while the black (dashed) line uses the tree-level relation. Shaded regions in 
both figures correspond to a stau mass below the LEP bound of 90 GeV. 

the associated value of y T to be the one read from the line labeled "Loop", giving « 0.585. 
In addition, for this particular choice of m L3 ,m E3 , /i, A T and M 2 , this figure also provides 



a translation between tan/3 and tan/3 re ff, as defined in Eq. (12). For instance, as stated 



above, the value of the radiatively corrected y T associated with tan /3 = 50 is ~ 0.585. The 
tan/3 TC fj corresponding to this y T can then be read-off to be ~ 57.5 using the line labeled 
"Tree" (which is simply the relationship defined in Eq. (|7|). 

From the above analysis, it is clear that including the effect of A T allows for a larger 
tau Yukawa coupling, y T , which increases the values of the stau mixing parameter, fiy T , 
allowed by metastability constraints. Such large stau mixing effects, in turn, allow for a 
larger enhancement of the rate of the Higgs decay into diphotons. 



B. Higgs Couplings to rr and bb 

As mentioned briefly before, smaller values of rriA and sizable values of A T induce an 
additional new physics effect in the diphoton event rate associated with a reduced Higgs 
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total decay width. In the MSSM, the lightest CP-even Higgs is a linear combination of h u 
and hd. The mixing angle, a, is governed by the off-diagonal element of the CP-even Higgs 
mass matrix, 



M 2 



H 



m\ sin 2 /? + Mf cos 2 f3 - (m 2 A + M§) sin /3 cos /3 + Loop 12 
(m 2 A + M|) sin cos + Loop 12 cos 2 /3 + M| sin 2 (3 + Loop 22 , 



(15) 



where [2 |20] 



Loop 



and 
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16tt V sin 2 /3 M 2 USY 



4 2 



M 2 

iw SUSY 



. 2 jj?Ab y$v 2 . 2 



16tt 2 



iw SUSY 



487T 2 



Mi 



sin(2a) 



2(A<|r)i2 



, (16) 



(17) 



^Tr[M 2 H ] 2 -4det[M 2 H ] 
At tree level, h u couples only to the up-type fermions and hd to the down-type fermions 
(leptons and down-type quarks). Since the lightest Higgs is given by the combination 



h = —hd sin a + h u cos a , 

its tree-level coupling to down-type fermions are then given by 

sin a rrid 



9hdd 



cos (3 v 



Furthermore, Eq. (19) gets corrected at one-loop 

sin a rrid 1 



9hdd 



1 



(18) 



(19) 



(20) 



cos (3 v (1 + A d ) \ tan a tan (3 y 

In the decoupling regime where rriA is large, we have sin a — > — cos/3 and cos a — > sin/?, 
so that the lightest CP-even Higgs couplings to fermions approach their SM values, but 
significant departures from these values may be obtained for smaller values of tua- In the 



absence of loop corrections to the Higgs mass matrix elements, Eq. (16), the down- type 



fermion couplings to the Higgs tend to be enhanced with respect to the SM values for 
moderate or small values of ttia [21] • However, for a non-zero and positive fiA T , the loop- 
corrections may lead to a relevant suppression of the off-diagonal term in the Higgs mass 



matrix in Eq. (15) and consequently a reduction of | sina|. As a result, the bottom and r 



couplings of the lightest CP-even Higgs boson may be also suppressed. The outcome is that 
the Higgs total width may be reduced and the branching fractions into gauge boson pairs, 
including diphotons, will be enhanced. 



11 



Note that for large values of A T , the suppression in the hrr coupling is larger than the one 
in the hbb coupling. This follows from the fact that the value of sin a is small and negative 
and sin a/ cos (3 ~ tan a tan (3. For \A d \ < 1, a reduction of the Higgs decay into down-type 
fermions may only be obtained for | sin a/ cos/3| < 1. This then implies that the coefficient 



of Ad in the numerator of Eq. (20) must be positive and larger than one. Therefore, a 
negative (positive) would decrease (increase) the ghdd coupling compared to the case 
A d = 0. Since in our scenario A b is positive and A T is negative, a sizable suppression of the 
r coupling of the Higgs may be induced for large values of tan f3 and positive A T , while the 
Higgs coupling to bottom quarks remains closer to the SM value. 



C. Results 



In the following we will consider tan/3 re g as an input parameter to study the interplay 
between the vacuum stability constraint and the possible enhancement of the Higgs to 
diphoton partial width. 

In Fig. [2] we show the bound on \x tan /9 Te ff from vacuum stability constraints as a function 
of the lightest stau mass, for A T = 0, = 2 TeV, tan/3 re ff = 70 and 6h — 1- In the left 
panel, we show the "absolute stability" bound obtained by imposing that the electroweak 
minimum is the global minimum of the theory; in the right panel we show the "metastability" 
bound obtained by imposing that the electroweak minimum is only a local minimum but 
with a life time longer than the age of the Universe. The corresponding bounds from Ref. [13] 
are shown by the black dashed lines. 

We compute the metastability bound using the following two methods: 

• The numerical package CosmoTransitions [27J, which uses a path deformation tech- 
nique to compute the bounce solution for a multi-dimensional scalar potential. We 
refer the reader to Ref. [28] for details. The result corresponding to the scalar potential 
in Eq. is shown by the blue solid line. 

• A numerical procedure which computes the bounce action by reducing the problem to 
a one dimensional one: at large values of /ztan/3, we define a canonically normalized 
field, $, which connects the charge breaking minimum and the electroweak breaking 
minimum. We compute the scalar potential of the new field, Vi($), using the potential 
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FIG. 2: Right panel: Metastability bound on /itan/3 re g- as a function of m^, having fixed A T = 0, 
m,A = 2 TeV and tan/3 e g- = 70. The black dashed line is the bound obtained in Ref. }13f . the blue 
solid line is the bound obtained by CosmoTransitions and the red solid line is the bound obtained 
solving the one dimensional equations of motion (see text). Left panel: Analytic absolute stability 
bound on /xtan/3 re ff as a function ofm^, for the same set of supersymmetry parameters as in the 
right panel, and the comparison with the bound from Ref. |13j . 

in Eq. Q. In this one-dimensional setup, the bounce solution is then calculated using 
the conventional "over-shoot/under-shoot" method [29]. The red solid line shows this 
result. 

We emphasize that, in order to compare with Ref. |13j . the results of the above two methods, 
shown in the right panel of Fig. |2j are evaluated using the improved tree- level scalar potential 
m Eq. Q. One sees that the metastability bounds obtained by the two methods we use are 
in excellent agreement with each other, however, they are significantly less stringent than 
the bound from Ref. [13] 4 . While the metastability constraint involves finding the minimal 
path connecting the two minima and computing the resulting bounce action numerically, the 
absolute stability bound is unambiguous. Thus in the left panel of Fig. [2] we also compare 

4 To confirm our results, we used a third method based on Ref. [30], which approximates the one dimensional 
potential described above with a triangle. The resulting bound is again in good agreement with the other 
two methods we used. 
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FIG. 3: Left panel: Metastability bound on //tan/3 re g- as a function ofm^, formi z = niE 3 , A T = 0, 
m A = 2 TeF, Mi = 55 GeV, M 2 = 400 GeV and M 3 = 1200 GeV. Right panel: Enhancement 
in the diphoton partial width with respect to the SM expectation, as allowed by the metastability 
condition, as a function of the lightest stau mass and for the same supersymmetry parameters as 
in the left panel. 

the absolute stability bound, computed analytically, with the one obtained in Ref. [13] 5 , and 
again find differences similar to those obtained from the metastability bound comparison 6 . 

Our main goal is to study the tan (3 dependence of the vacuum metastability bound on 
/xtan (3 Te ft- CosmoTransitions is a public code with the capability of handling the full scalar 
potential, encompassing both the up- and the down-type Higgses, as well as the one-loop 
effective potential. Therefore, in the following we choose to present our results based on 
the outcome of CosmoTransitions using the full one loop effective potential, imposing stop 
mass parameters consistent with a Higgs mass of about 125 GeV. It turns out that using 
the one-loop effective potential instead of the improved tree-level potential only results in 
a difference of a few percent on the bound on fi tan (3, which explains the small differences 
between Figs. [2] and |3j 

In the left panel of Fig. [3j we present the bound on \i tan (3 Te s as a function of the lightest 

5 We extrapolated the results presented in Fig. 4 of Ref. [T3] . 

6 See footnote [lj 
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FIG. 4: Ze/i panel: Metastability bound on /itan/3 Te g as a function of m^, /or mj, 3 = me 3; 
Ar = 1 TeF, = 1 TeF, Mi = 55 GW ; M 2 = 400 GeV and M 3 = 1200 GeV Right panel- 
Enhancement in the diphoton partial width with respect to the SM expectation, as allowed by the 
metastability condition, as a function of the lightest stau mass and for the same supersymmetry 
parameters as in the left panel. 

stau mass for different choices of tan/3 re fj. We set m Lz = m E3 , Mi = 55 GeV, M 2 = 400 
GeV and M 3 = 1200 GeV. This first plot shows the results in the decoupling limit with no 
CP-even Higgs mixing from the stau sector: m A = 2 TeV and A T = 0. One sees that the 
bound becomes weaker as tan/3 re ff grows. As explained in Section [IT], this is because for a 
fixed value of fj,ta.nf3 Te s the stabilizing quartic term It^-trI 2 increases with tan/3 2 efr . In the 
right panel of Fig. [3] we show the allowed enhancement in the diphoton partial width for the 
same set of parameters. One sees that an enhancement of up to 50% may be obtained for 
tan/Jreff < 100 and a stau mass of 90 GeV, the LEP limit. Larger enhancements may be 
achieved for even larger value of tan (3 Te s 7 . 

7 The value of Mi we chose is of the right order to generate a proper Dark Matter relic density for values of 
the lightest stau mass of about 100 GeV [5]. For such values of Mi (2Ml < m/j), a small invisible width 
is generated, which becomes more significant for larger values of tan j3 and A T . This enhancement in the 
invisible width is due to an increase of the lightest neutralino Higgsino component with decreasing values 
of /i. The total width is then enhanced by at most a few percent in the region of parameters under study, 
and therefore leads to a reduction of all visible branching ratios by a similar amount. 
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In Fig. [4] we show similar plots with A T = rriA = 1 TeV. We again set = rriE 3 , 
Mi = 55 GeV, M 2 = 400 GeV and M 3 = 1200 GeV. Comparing the left panels of Fig. [| 
and Fig. [|J we note that the bound on /itan/3 eff is about 20% more stringent for the lower 
value of rriA and larger value of A T . As discussed before, this is due to the destabilizing 
effect of the A T trilinear coupling in Eq. ([9]), which can generate new charge breaking vacua 
at relatively large and negative values of the field hd- 

As can be seen from comparing the right panels of Fig. [3] and Fig. [4], the effect on the 
Higgs diphoton decay rate of the more stringent bound on /ztan/3 e ff, in the case of larger 
A T and smaller tua, may be partially compensated by the suppression of the bb width as 



discussed below Eq. (20): an enhancement of the diphoton rate by a factor of ~ 40% may 
still be obtained for tan /3 TC ff < 100 and a lightest stau mass at around 90 GeV. 

In order to compute the branching ratio of the Higgs decay into bottom-quarks and tau- 
leptons we used the latest version of the public program CPsuperH. This program computes 
the Higgs spectrum and decay rates, as well as the sparticle spectrum, including the A r ^ 
effects, and therefore provides a consistent framework for computing the Higgs decay widths 
as a function of the lightest stau mass within the effective theory approximation described 
in this work 8 . The suppression of the Higgs to tt and bb rates is shown in Fig. [H] It is clear 
from this figure that for these values of A T and tua, a relevant suppression of the tt rate is 
only possible for very large values of tan/3 rcff . 

Larger values of A T may lead to larger stau contributions to the off-diagonal term of the 



CP-even Higgs mass matrix element, Eq. (16). However, larger values of A T also induce 



a stronger metastability bound on /itan/3, which in turn implies that the effect of A T on 
the Higgs mixing is reduced. We checked that the combination of all these effects is such 
that varying the value of A T within a few hundred GeV leads to only small changes on the 
results presented in Fig. |5j For increasing values of A T , we found that the suppression of 
\i tan j3 was such that overall there was smaller allowed suppressions of the Higgs to bb and 
Higgs to tt decay branching ratios. Similar comments apply to variations of rriA, although 



A quantitative comparison of the results obtained from CPsuperH with the ones obtained from FeynHiggs 
[5Tj is technically difficult, since FeynHiggs does not include A r effects, but does include effects as 
well as full one-loop corrections to the /irr coupling. Using FeynHiggs with taji/3 Te g and A Te ff as input 
yields a similar stau spectrum as the one from CPsuperH with the corresponding tan /? and A T , but the 
bottom Higgs coupling is artificially affected by this change. 
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FIG. 5: Branching ratio of the decay of the lightest CP-even Higgs into (Left): two r leptons, and 
(Right): two b squarks as a function of the lightest stau mass, normalized to their SM values, for 
m L3 = m Ez) A T = 1 TeV, m A = 1 TeV, M x = 55 GeV, M 2 = 400 GeV and M 3 = 1200 GeV 

in this case it is smaller values of m A that may lead to larger CP-even Higgs mixing effects, 
which are limited by stronger bounds on //tan/3. Therefore, Fig [5] is representative of a 
more general case and, quite generically, very large values of tan /3 Te g are necessary in order 
for stau effects to modify the Higgs to bb and Higgs to rr decay rates in a relevant way. 



III. CONSTRAINTS ON LARGE TAN (3 

The large tan/3 re a needed to enhance the diphoton width (and suppress the rr/bb cou- 
plings) leads to large values of y T and yb at the weak scale. One may be concerned that 
a large tau/bottom Yukawa coupling may become non-perturbative at some intermediate 
energy scale below the GUT scale, necessitating new physics beyond the MSSM at or below 
the Landau pole energy. 

To get a better sense of the values of tan /3 and tan /3 Te ff in the region of interest, in Fig. [6] 
we plot A r and y T as a function of tan /3, similar to Fig. [TJ but this time taking into account 
the metastability constraints on the // parameter. Although in the figure we have chosen 
specific values of the soft supersymmetry breaking slepton masses, we verified that A T is 
not very sensitive to the resulting stau masses. Therefore, the right panel of Fig. [6] can be 
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FIG. 6: Left panel: A T versus tan /3 for m L3 = m E3 = 240 GeV, Mi = 55 GeV, M 2 = 400 GeF, 
M3 = 1200 GeF and A T = 1 TeF. T/ie lightest stau mass varies in the range 92-103 GeV, with 
the lighter masses associated with the larger values o/tan/3. The values of fi are taken to be those 
associated with the metastability limit. Right panel: y T versus tan/3 for the same set of parameters. 



The red (solid) line uses the one-loop relation (Eq. (10)), while the black (dashed) line uses the 
tree-level relation (Eq. Pfy)- 

used to approximately translate between values of tan /3 re g and tan (3 for the entire range of 
parameters of interest in this work. 

The running of the r and bottom Yukawa couplings depends strongly on their values at 
the weak scale, which are determined by tan/3 re ff and the analogous tan/3& e ff, respectively. 
In the region of parameters we are studying, namely soft breaking squark and gluino masses 
of order 1 TeV and stop mixing A t of about 1.5 TeV, one obtains A& ~ 50%. This leads 
to values of the bottom Yukawa coupling at the TeV scale of order 0.6 - 1, rather than 
the values of order 0.9 - 1.5 that would have been obtained had we used the tree-level 
relationship for values of tan/3 ~ 60 - 100. In addition, the effect of the QCD coupling 
on the evolution of the quark Yukawa couplings decreases their values for increasing energy 
scales. Thus, rather than from the bottom Yukawa coupling, the strongest perturbativity 
constraint comes from the running of the r Yukawa coupling, due to its large values at the 
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FIG. 7: Two-loop RG evolution of the r (solid lines) and b (dashed lines) Yukawa couplings, for 
different values o/tan/? Te g. We terminate the running of the bottom Yukawa at the same scale at 
which the tau Yukawa for the corresponding tan/3 re g becomes non-perturbative. 

weak scale induced by a negative A T . 

Fig. [7] shows the two-loop renormalization group (RG) evolution (32] of y T and as 
a function of the RG scale Q, using the weak scale Yukawa couplings obtained from in- 
cluding A T and A;, effects. We see that perturbative consistency up to the GUT scale, 
Log 10 [Q/GeV] = 16, may be obtained for tan/3 re fr < 90, which corresponds approximately 
to tan/3 < 70 according to Fig. [6} Larger values of tan/3 Te g demand an ultraviolet comple- 
tion at scales below the GUT scale due to the appearance of a Landau pole. Comparing 
these results with the ones presented in Fig. [5j we observe that, in the light stau scenario, 
a suppression of the Higgs decay branching ratio into r leptons larger than ~ 5% would 
require an ultraviolet completion of the MSSM at scales below the GUT scale. 

Another constraint on the very large tan regime comes from flavor physics. In the 
Minimal Flavor Violation hypothesis [34J, the most important flavor observables receiving 
tan/3 enhanced new physics contributions are B u — > tu, b — > 57 and B s — > fifi. These 
observables could in principle give a stringent bound on the value of //tan/3. It has been 
shown recently that, assuming A t > and stop masses of about 1 TeV, the b — > S7 branching 
ratio is enhanced, and the B s — > fifi tends to be smaller than the SM expectation [36J. 
However, although consistency with the observed B s — > fi^i value is more easily achieved for 
fiA T > 0, quite generally the bounds coming from b — > 57 and B s — )■ /i/i depend strongly 
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on the splitting between the third generation and the rest of the squark masses, an effect 
not directly related to the Higgs decay rate into diphotons ( see, for example, Ref. |35J, for 
a recent discussion). A more robust constraint is instead represented by B u — > ru: charged 
Higgs contributions to this decay rate arise at the tree level and are generically large for 
large values of the tau and bottom Yukawa couplings. It has been shown, however, that, 
choosing a charged Higgs mass at around 1 TeV would satisfy the bound from B u — > tu, for 
the entire range of values for y T and considered in this paper [37]. 

Finally, we note that MSSM heavy Higgs bosons at around 1 TeV are starting to be 
probed by direct searches through their decays into rr and bb [33]. A naive extrapolation 
of the CMS results presented in Ref. [33] would indicate that values of rriA — 1 TeV may be 
ruled out a the 95% C.L. for values of tan/3 > 70. However, for such large values of tan/3 
and A T , the width of these Higgs bosons is very large. In addition, these heavy Higgs bosons 
have relevant decays into stau pairs, which may suppress the branching fraction into the 
tau leptons [39J. Therefore, a naive extrapolation may not be valid and a detailed analysis 
is necessary in order to determine the constraints on the CP-odd Higgs boson mass from 
direct searches. 

Other phenomenological constraints, like precision electroweak observables, the anoma- 
lous magnetic moment and the dark matter relic density in the light stau scenario were 
discussed in Refs. [HI [38] with positive conclusions. These previous studies, however, did not 
include the effects of A r , which are important in this region of parameter space. 

IV. CONCLUSION 

In this work we studied the vacuum stability constraint on the light stau scenario in 
the MSSM, which can produce a significant enhancement in the Higgs-to-diphoton decay 
width. We improved upon earlier studies of the metastability condition by analyzing in 
detail the tan/3 dependence and also including the effects of a non-zero A T , which corrects 
the relationship between the tau lepton mass and the tree-level tau Yukawa coupling. We find 
that an enhancement of the order of 50% is consistent with the requirements of metastability 
of the electroweak vacuum and the perturbativity of the MSSM up to the GUT scale. 

In the region of parameter space we are interested in, A T is negative and of the order of 
10 - 25%. Hence, the tree-level tau Yukawa, given the measured tau lepton mass, becomes 
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larger than when neglecting A T . Since the stabilizing quartic term of the stau-Higgs scalar 
potential, \f L fji\ 2 , is proportional to the square of the tau Yukawa coupling, y 2 , a larger y T 
relaxes the bound coming from vacuum stability. Such an effect becomes more important 
for larger values of tan (3, and hence it is relevant to take into account the tan (3 dependence 
on the metastability constraints on //tan/3. 

In addition, we also studied the impact of having a non-zero A T on both the vacuum 
stability constraint and the Higgs to diphoton decay width. We found that positive values 
of \xA T may significantly impact the metastability condition, lowering the bound on // tan (3 
by ~ 20% for A T = m A = 1 TeV. On the other hand, we also showed that A T impacts 
the CP-even Higgs mixing and thus can decrease the Higgs decay into taus and bottoms 
depending on the value of uia- In particular, for sizable A T and fiA T > 0, the suppression of 
the Higgs to diphoton width due to lower allowed values of // tan (3 is partially compensated 
by the decrease in the total Higgs width, leading to a reduction of no more than ~ 10% in 
the Higgs decay rate into diphotons for A T = m A = 1 TeV with respect to the A T = case. 

We then studied the two-loop RG running of the tau and bottom Yukawa couplings and 
found that an enhancement in the diphoton width of the order of 50% is consistent with 
perturbative values of the Yukawa couplings up to scales of the order of the GUT scale. 
On the contrary, a significant suppression of the width of the Higgs decay into tau leptons, 
larger than ~ 10%, requires an ultraviolet completion of the MSSM at scales below the GUT 
scale. 

Finally, though more data is necessary to determine if the enhancement suggested by 
current measurements is real or a product of statistical fluctuations, we have shown that 
the precise Higgs decay rate to diphotons in the MSSM is intimately connected to the fate 
of the Universe. In conclusion therefore, a precise measurement of the Higgs coupling to 
diphotons has far reaching implications, and should be a high priority at the LHC and at 
any future Higgs factory. 
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APPENDIX 



In this appendix we provide, for convenience, the approximate analytic expressions for 
A T and Aj,, the complete expressions given in Ref. [20]. We first define the loop function 

^ _ a 2 b 2 \og(a 2 /b 2 ) + b 2 c 2 \og(b 2 /c 2 ) + c 2 a 2 log(c 2 /a 2 ) 
Then in the region tan ^ > 1, one can write 



3/i tan 8 „ , , ,, , , , a tan B , , , r . 
A T ^ - P QO 2 9 2 2 M 2 I(m Pr ,M 2 , fj.) - y~9\ M x /(m fl , m f2 , M 1 ) , (22 

3Z7T Z lD7T 

/xtan (3 



^M- g I(m- bl ,m- b2 ,M d ) + ^A t I(m {l ,m i2 , fi) 



(23) 



2vr 2 

In the above Mu2) is the mass parameter for the Bino (Wino), the stau masses, 

Mg is the gluino mass, mi 2 are the sbottom masses, A t is the trilinear soft-breaking term 
in the stop sector, and 2 are the stop masses. Moreover, g x is the gauge coupling for 
U(1)y, gi is the gauge coupling for SU(2) L , g 3 is the gauge coupling for SU(3) C , and y t is 
the top Yukawa coupling. 
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